Elect romagnon excitations in modulated multiferroics 
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The phenomenological theory of ferroelectricity in spiral magnets presented in [M. Mostovoy, 
Phys. Rev. Lett. 96, 067601 (2006)] is generalized to describe consistently states with both 
uniform and modulated-in-space ferroelectric polarizations. A key point in this description is the 
symmetric part of the magnetoelectric coupling since, although being irrelevant for the uniform 
component, it plays an essential role for the non-uniform part of the polarization. We illustrate 
this importance in generic examples of modulated magnetic systems: longitudinal and transverse 
spin-density wave states and planar cycloidal phase. We show that even in the cases with no uniform 
ferroelectricity induced, polarization correlation functions follow to the soft magnetic behavior of the 
system due to the magnetoelectric effect. Our results can be easily generalized for more complicated 
types of magnetic ordering, and the applications may concern various natural and artificial systems 
in condensed matter physics (e.g., magnon properties could be extracted from dynamic dielectric 
response measurements) . 

PACS numbers: 75.30.Ds, 75.50.-y, 75.25.-l-z, 62.12.Br 



The existence of magnetic materials that are also fer- 
roelectrics P, 0, Q upsurge recurrent interest to study 
and to rationalize theoretically various models manifest- 
ing this possibility (see e.g. 0, 0, H, 01 )■ The recent 
discovery of ferroelectricity in a family of the so-called 
frustrated magnets, see e.g. Q, is particularly instruc- 
tive. In this case, ferroelectricity is a collateral effect of 
the magnetic ordering that, in contrast to earlier multi- 
ferroics, involves modulated structures (ol. \ld^. The cur- 
rent interest in this new class of ferroelectromagnets is 
twofold. On one hand the fundamental physics that re- 
sult from the interplay between the electric and magnetic 
properties in these systems is not exhausted [loj and, on 
the other hand, many potential applications arise from 
this interplay 



The appearance of a modulated magnetic structure 
is accompanied by a distribution of polarization due to 
the inhomogeneous magnetoelectric effect [12} • This uni- 
versal effect is similar to so-called flexoelectric effect in 
liquid crystals [l3|- Microscopic studies have focused 
on the Dzyaloshinskii-Moriya interaction as the main 
source of the inhomogeneous magnetoelectric effect 
and this idea seems also to be behind the phenomeno- 
logical Landau-like approaches carried out until now in 
relation with ferroelectricity of magnetic origin 
As a result, the main attention has been put to only an- 
tisymmetric parts of the effect (where the components 
of the magnetization enter in an antisymmetric combi- 
nation). However, in a general case, the inhomogeneous 
magnetoelectric effect has also a non- vanishing symmet- 
ric part (see below). This is indeed the natural possibility 



from the phenomenological point of view 
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what- 



ever the microscopic origin of this symmetric coupling. 

In this paper we show that this symmetric part of the 
inhomogeneous magnetoelectric effect is a key ingredi- 



ent in the full description of materials exhibiting mod- 
ulated magnetic structures. When it comes to the uni- 
form polarization, indeed, only the antisymmetric part of 
the magnetoelectric effect plays a role in the agreement 
with previous theoretical publications (see e.g., P. Isl. ITst . 
and more details below). However, to describe non- 
uniform, i.e. modulated in space, polarizations the mag- 
netoelectric effect has to be considered in its full extent. 
This is precisely the case when addressing the so-called 
clcctromagnons, i.e., the normal modes characteristic of 
ferroelectromagnets that involve both polar modes and 
magnons in a hybridized way 0, • 

Electromagnons were theoretically predicted long time 
ago considering uniformly ordered materials This 
consideration has been recently extended to the case of 
an antiferromagnet that becomes ferroelectric and then, 
through the inhomogeneous magnetoelectric effect, this 
induces an incommensurate magnetic structure [it} . The 
first experimental evidences about electromagnons have 
been reported for materials in which, showing no instabil- 
ity towards ferroelectricity in the absence of magnetism, 
it is the appearance of a modulated magnetic structure 
what make them ferroelectrics 18, [l^. Unfortunately, 
all theoretical studies performed so far for this latter case 
have been restricted to oiily antisymmetric parts of the 
magnetoelectric coupling 14, 2^. We show below that, 
contrary to what is implicitly assumed in these works, 
both symmetric and antisymmetric parts of the inhomo- 
geneous magnetoelectric effect contribute nontrivially to 
the interplay between the fluctuations of magnetization 
and polarization in multiferroics. 

First of all, let us reconsider the uniform polarization 
that may arise due to the inhomogeneous magnetoelectric 
effect. To simplify our presentation and ease the algebra 
we neglect crystalline field anisotropy (a more realistic 
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treatment does not affect our qualitative results) . In this 
case, tlic magnctoelectric coupling term reads as 

Fem = -fiP ■ M(V • M) - /2P • [M X (V X M)], (1) 

where P is the polarization and M is the magnetiza- 
tion, and /i and /2 are, in general, two different con- 
stants. The general form of this term is — fki,ij PkMidiMj 
[l^ . where the magnetoelectric tensor reduces to fki,ij = 
fiSkiSij + f2{5kAj - SkjSii) in the isotropic case. Let 
us split the polarization into uniform P and non-uniform 
(modulated) P' contributions: P = P + P'. From the 
magnetoelectric term, we then can extract the following 
contribution to the free energy of the system [2l|: 

-fki,ijPk JdrMidiMj = ^ J dr{fkiji ~ fki,ij)MidiMj . 

Thus, we see, the uniform polarization indeed couples 
to the magnetization only through the antisymmetric 
part of the magnetoelectric tensor. In the isotropic 
case that is fki,ij - fkiji = (/i + f2){Ski6ij - SkjSu). 
So, as regards the uniform polarization, the relevant 
part of the magnetoelectric term ([1]) can be written as 
-^(/i + /2)P • [M(V • M) - (M • V)M], which coincides 
with the expression given in [5|. However, it is important 
to keep in mind that there are actually two independent 
contributions, /i and /2, to the magnetoelectric coupling 
term. 

Let us now turn our attention to the magnetic fea- 
tures of the modulated magnets discussed in e.g. [lo| . 
As we have mentioned, in these systems ferroelectricity 
is induced by magnetism and does not appears sponta- 
neously by itself. In addition, the magnetoelectric effect 
is typically weak in comparison with pure magnetic con- 
tributions into the free energy. Therefore, in the first 
approximation, when determining the magnetic structure 
and dynamics it is possible to consider the magnetization 
separately. 

Most of these systems first undergo a transition from 
the magnetically disordered phase to a modulated struc- 
ture similar to a longitudinal spin-density- wave (LSDW) : 

M = (Micosga;,0,0), (2) 

and then to structure similar to a planar cycloid (PC): 

M = (Ml cos Qx, M2 sin Qx, 0). (3) 

The dependence on P of Landau free energy density can 
be taken as \x^^P^ + Fem^ where x is the dielectric 
susceptibility in the absence of magnetism [2^ and Fem 
accounts for the magnetoelectric coupling. It can be eas- 
ily seen that the inhomogeneous magnetoelectric effect 
([T|) implies the appearance of ferroelectricity as a result 
of the PC structure, but not of the LSDW [5|. 



A simple way to reproduce this sequence of transitions 
is by considering Landau free energy density in the form 

Fm = ^A/' + + i ^ M,Um,, (4) 

i 

where the differential operators Li describe the 
anisotropic softening that gives rise to the modulated 
structures. For small wavevectors one can approximate 
Li — cV^. The free energy ([4]) then is reduced to 
a familiar form (see e.g. [23[). Near the wavevector 
Q = (Q,0,0) of the modulated structures a + Li ~ 
cti + CxV'^ + cV5_. Eq. dl]) is then a natural general- 
ization of the free energy considered in [l^ that pro- 
vides the unified description of the spectrum at the rel- 
evant wavevectors (i.e., q = and ±Q). If a,b,c > 0, 
the system shows no instability with respect to uniform 
magnetizations. Furthermore, if ct j, < ay < a^, in the 
spirit of arguments presented in one can see that 

the free energy is minimized by the LSDW modulation 
in the parameter range < 0, ay > 0, while the PC 
configuration occurs when < and 3ax ^ ay < ax/S 
(with az > in both cases). 

At this point, it is worth making the following com- 
ment on the label "improper ferroelectrics" that is fre- 
quently put to the systems of our interest (see e.g. \ld^). 
Within above scenario My or, more precisely, its ±Q 
Fourier components, can be seen as the order parameter 
of the LSDW-paraelectric ^ PC-ferroelectric transition. 
In the LSDW-paraelectric state Eq. ([T|) has a term cou- 
pling My and Py linearly so, in view of this linearity, 
the system should be labeled as proper rather than im- 
proper ferroelectric . This is in fact manifested in the 
corresponding anomalies Q . 

As regards the dynamics, we are interested in the low- 
frequency excitations of the system. In the magnetically 
ordered phases, not too close to the transition points, 
these excitations are associated with variations of the 
magnetization in which its modulus is conserved. Thus 
the excitations are described by the Landau-Lifshitz 
equation psj : 

M = 7M X Heft, (5) 

where 7 is a constant and Hcff = (L^M^, LyAIy, L^M^). 

The dispersion relations and correlation functions of 
our interest are obtained by linearizing this equation 
about the corresponding configurations of equilibrium. 
Then it is convenient to consider first of all the LSDW 
structure ([2]) , and then a virtual transverse spin-density- 
wave (TSDW) modulation: 

M = (0,M2sinQa:,0). (6) 

This would be the structure obtained from (jl]) if, con- 
trary to what we assumed for the LSDW, ay becomes 
negative with a^; > 0. This structure TSDW is of in- 
terest in its own right and besides, it turns out that the 
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low-lying normal modes of the PC case can be found as 
the superposition of the normal modes associated with 
LSDW and TSDW structures. This notably simplifies 
the corresponding calculations. 

Skipping a large amount of tedious (although straight- 
forward and simple) algebra, we come to the following 
results. They are obtained under the same assumptions 
as those for the static problem. That is, assuming that 
the anharmonicity is weak and that the spatial dispersion 



such that the generalized stiffness presents well defined 
minima only at the wavevectors q = 0, ±Q. This allows 
us to decouple the modes q and q ± Q from the rest [i^ . 
To our purposes, higher-order satellites can be neglected. 

As it could be expected, the low-lying modes of the 
LSDW represent independent oscillations of the magne- 
tization along the y and z axes. The corresponding dis- 
persion relations are such that 
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Here mx = 7-M1/2 and the quantities on are defined such 
that LiC^^"^^^ = one^^"^^^. It is worth noticing that the 
closer the transition is, the smaller is the gap obtained at 
small wavevectors. Not only because of the smallness of 
TOj,, what influences trivially on the whole spectrum, but 
also because Ux 0. In fact, very close to the transition 
point LOy^z{q Q) ^ q. The gap obtained close to the 
wavevector of the modulation, on the contrary, is not so 
sensitive to the smallness of ax ■ 

Similar expressions are found for the TSDW structure 
([6|) . In this case, the low- lying modes are associated with 
oscillations of the magnetization along the x and z axes. 
Their dispersion relations can be obtained from ^ by 
replacing y x in nix, cti and S^. 

As we have mentioned, these results for the LSDW 
and TSDW structures merge in the PC case. Three low- 
lying normal modes are found such that the modulus 
of the magnetization is conserved in this latter struc- 
ture. The modes associated with oscillations along the 
X and y axes coincide with that obtained previously for 
the TSDW and LSDW modulations respectively. Con- 
sequently, they have the same dispersion relations. The 
mode representing oscillations along the z axis, however, 
is composed by the modes associated with the same type 
of oscillation in the TSDW and LSDW cases. In conse- 
quence, its dispersion relation can be expressed such that 



(q) 



■'^.LSDW 



(q) 



•'z.TSDW 



(q)- 



Let us now turn our attention to the fluctuations. 
Since the softness of the system is in its magnetic part, 
fluctuations of the magnetization will play the main 
role. Once Eq. ([5]) is linearized, the corresponding cor- 
relation functions: 



Af,,(q) = {SM,iq,Lu)SM,i^q,-Lu)), 



(8) 



can be found with the aid of the fluctuation-dissipation 
theorem [2^. Within a first approximation, the polar- 
ization will follow these fiuctuations in accordance with 
the magnetoelectric coupling ([T|). This results in fluc- 



tuations of the polarization SP of magnetic origin, i.e., 
electromagnons, whose correlation functions: 



Py(q) = (<5P,(q,c^)JP,(-q, -u)), 



(9) 



can be expressed in terms of the previous quantities ([8]). 
The hybridization between magnons and polar modes 
reflects also in the form of cross correlation functions 
(SPiSMj). Within our approximation, these cross cor- 
relations can also be expressed in terms of the quantities 
([8|). But since they essentially reveal the same informa- 
tion (see e.g. [13), we omit them in the following. 

As regards magnetic correlation functions, the "non- 
diagonal" components are zero for the structures of our 
interest (My = if i 7^ j). In addition, Mxx = in the 
LSDW case while Myy = in the TSDW one. Close to 
the relevant wavevectors, the non-vanishing correlation 



functions are such that Afjj^(q) cx t{[u 



T '^uj'^}, where t is the (phenomenological) relaxation 
time of the corresponding oscillations. 

The fluctuations of the polarization of magnetic origin 
can be described in terms of the quantities Ma. These 
fluctuations are such that the soft magnetic behavior at 
q « ±Q reflects in dipole polarization excitations with 
small wavevectors. And reversely, the soft magnetic re- 
sponse at (7 — > reflects in polarization excitations with 
q w ±Q, as can be seen in further formulas. Note also 
that the magnetoelectric coefficients /i and /2 enter sep- 
arately in the following results. That is, not only through 
the combination /i -I- /2 that determines the antisymmet- 
ric part of the magnetoelectric effect and, consequently, 
the space average of the polarization of the system. (The 
q-dependence associated with x 123 is implicit in these 
formulas.) 
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In the LSDW case we have 

PccA^) - f!{qy [Myyiq - Q) + Myy{q + Q)] 

+ ql [M,,(q - Q) + RL^iq + Q)] }, 
Pvvin) - 9U{q)Myyiq - Q) + glM)Myy{q + Q), 
P.ziq) ~ 9\{q)M,,{q - Q) + gl^{(i)M,,{(i + Q), 
Pxy{q) ^ fiqy [9-2iq)Myy{(i - Q) - g+2{q)Myy{(i + Q)] , 
Px^iq) ~ fiQy [.g-2(q)Af..(q - Q) - .g+2(q)i\/..(q + Q)] , 
PyM = 0. 

where .g±2(q) = (/i + f2)Q ± /2'Z2;- In the case of a pure 
symmetric couphng (/i = —f2), with no Dzyaloshinksii- 
Moriya-hke contributions, all these correlation functions 
turn out to be oc q^. In consequence, the eventual im- 
portance of fluctuations of the polarization with small 
wavevectors is conditioned to the existence of a non- 
vanishing antisymmetric part in the magnetoclcctric cou- 
pling. In any case, fluctuations with wavevectors q « ±Q 
arc important irrespective to the above vanishing (they 
arc entirely due to /i 7^ 0). These "magnetic" softenings 
of optical phonons in the LSDW paraelectric state can 
be used to measure the corresponding coefficients /i and 
/2. This soft behavior in the absence of ferroelectricity 
has in fact been observed in [l9| . 

Similarly, for the TSDW structure we find that 

^^..(q) - fWy [M^q. - Q) + A4x(q + Q)] , 
Pyyi^) ~ .9-i(q)M,,(q - Q) + .g^ ^(q)A4,(q + Q) 

+ flql [M,M - Q) + + Q)] , 

Pzz{q) - fhl [Mzziq - Q) + + Q)] , 

Pxy{q) - /2% [.g-i(q)A4,(q - Q) - ,g+i(q)A4,(q + Q)] , 

P..(q) = 0, 

PyM) - -hf2qyqz [M,,(q - Q) + + Q)] , 

where g±i(q) = (/i + f2)Q ± fiqx- As in the previ- 
ous case, the relevance of fluctuations of the polarization 
with small wavevectors is made conditional on the anti- 
symmetric part of the magnetoclcctric coupling while the 
relevance of fluctuations with wavevectros close to that 
of the modulated structure is not. 

Finally, armed with the knowledge of correlation func- 
tions for LSDW and TSDW structures, we are in the po- 
sition to calculate the polarization correlation functions 
in the PC case as Pl-^^"^ + P^j^^"^ + APij . Here the first 
two terms represent the correlation functions obtained in 
for the LSDW and TSDW modulations respectively, and 



the third one is 

APxx(q) ^ ql [/f n,(q, Q) + /|n^(q, Q)] , 

AF^^(q) ^ ql [/|n,(q, Q) + /fn^(q, Q)] , 
AP..(q) ^ fWz [n,(q, Q) + n,y(q, Q)] , 
AP,y(q) ^ /i/2feqJn,(q,Q) +n,(q,Q)], 
AP,,(q) ^ /i/2fe<7. [n,(q, Q) -t- n,(q, Q)] , 
AP^,(q) ^ hf2qyqz [n,(q, Q) + n,y(q, Q)] . 

where n,(q, Q) = mf [A/jj(q - Q) 4- Mu{q + Q)]. These 
additional contributions are relevant for the fluctuations 
with wavevectors close to ±Q. It is worth noticing that, 
within our approximations, there are no contributions 
oc miTO2 to these correlations functions. 

In summary, we have presented the key aspects of the 
phenomenological theory of ferroelectricity of magnetic 
origin that allow one to treat consistently, on the same 
footing, uniform and non-uniform polarizations. One of 
its key points is the accounting for all the possible contri- 
butions, symmetric and antisymmetric ones, to the inho- 
mogeneous magnetoclcctric effect. Considering the most 
representative modulated magnetic structures, we have 
employed this theory to study the corresponding fluc- 
tuations of magnetization and polarization. Our main 
findings are the following, i) The aforementioned fluctu- 
ations turn out to be interdependent as long as magnetic 
order appears, irrespective of whether this order induces 
an uniform polarization or not. This implies that ferro- 
electricity may not be strictly necessary when addressing 
applications based on such an interplay, ii) The softness 
of the system reveals in the form of large fluctuations at 
both, small wavevectors and wavevectors close to that 
of the modulated structure. At small wavevectors, the 
fluctuations of the polarization of magnetic origin can 
be associated with the antisymmetric part of the mag- 
netoclcctric effect (that eventually gives rise to the uni- 
form polarization of the system). At wavevectors close 
to that of the modulated structure, however, the situa- 
tion is different. In this case, the fluctuations are due 
to both symmetric and antisymmetric parts of the effect. 
Consequently, in view of the relativistic nature of the an- 
tisymmetric part, one can expect here a leading role of 
the symmetric contribution (which has its origin in the 
exchange interaction). 
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